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Photons can gravitationally scatter off a cosmic string loop and gain or lose energy. We consider the 
spectral distortion induced by cosmic string loops placed in an ambient thermal bath of photons. 
The fractional deviation from a thermal spectrum caused by cosmic strings is estimated to scale 
as (Gnm)^z^, where Gn is Newton's constant, fi is the string tension, and z is the cosmological 
redshift after which spectral distortions can survive. This effect is large enough to potentially be of 
observational interest. 



I. INTRODUCTION 

The spectrum of cosmic microwave background (CMB) 
photons has been measured very accurately and, at the 
present time, shows no deviations from a blackbody dis- 
tribution. These measurements provide tight constraints 
on any source of energy injection into the primordial 
medium at redshifts z < 10®. 

The relativistic motion of cosmic strings may con- 
tinually inject energy into the cosmic medium starting 
from very early epochs, and certainly much earlier than 
z — 10®. Depending on the particle physics interac- 
tions between the string and the cosmic medium, the 
energy deposited could be significant [1 . However, there 
is also energy deposition due to the gravitational in- 
teractions of photons and strings, and this process is 
completely generic. The gravitational scattering of pho- 
tons by strings can result in gravitational lensing of 
background structure with angular separation ^ AirGjsifi 
where Gn is Newton's constant and /i is the string's lineal 
energy density. However, the bending of light is an elastic 
process and cannot lead to spectral distortions. Photon 
energies can only be re-distributed by cosmic strings at 
order (G^^)^ and this is the effect we are interested in 
analyzing. 

Inelastic scattering by strings is similar to photon pair 
creation by the time-dependent metric of a cosmic string 
as discussed in [2]. Indeed the Feynman diagram is the 
same for both processes except for the reversal of the mo- 
mentum in an external photon leg. However, important 
technical differences arise due to the different kinematics 
that make the scattering calculation more difficult. To 
obtain the spectral distortions we also need to evaluate 
a Boltzmann integral which makes the analysis challeng- 
ing. 

In Sec.|TT]we introduce the formalism and evaluate the 
scattering amplitude for a general cosmic string loop. To 
obtain numerical results, we specialize to a simple loop 
trajectory, one that is expected to be close to the most 
generic loop [5], and is called a "degenerate, kinky loop". 
In Sec. |III| we numerically evaluate the photon scattering 
rate in phase space after averaging over all possible ori- 
entations and shapes of the degenerate kinky loop. The 
spectral distortion is shown in Fig. [2] and the result is 



used to estimate the fractional spectral distortions of the 
CMB in Sec. |IV[ In Appendix [A| we discuss some techni- 
cal details of the evaluation of the Boltzmann integrals. 



II. SETUP 

We will consider the process of photons scattering off a 
cosmic string taking place in a weakly curved spacetime, 



h 



Mpi = l/V327^G^ 



77^^ = diag[l, -1,-1,-1], 
The scattering process is governed by the action 



S = S. 




'NG 



(1) 

(2) 
(3) 



h^^Tl^^-^d x + J da j At /i^^T^'^'q) | , 



where is the action of the photon A^^ 
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(4) 



with F^^ = dfj^Ai^—d^Afj^; Sh is the quadratic action of the 
graviton hfj^^ in de Bonder gauge = (l/2)9y/i). 



(5) 



with h = rj^^'^h^y, and S'ng is the Nambu-Goto action for 
the cosmic string 



S'ng = -fJ- I do- / dt^/^ 



(6) 



Here, our cosmic string is of length L with tension /x, 
and [t, a) is the 4 coordinate vector of a given point 
(t, a) on its worldsheet. The determinant of the induced 
metric of the latter is 7; denoting 



we have 



= dtX^", (X^')' = d^Xf", (7) 

(8) 



j^x^x'^ - {^x-xy 
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The stress-energies of the photons and the string are, 
respectively, 

T^T) = -F'^'^F''^ + ^ri^'-F^pF'^^, (9) 



and 



^(NG) 



(10) 



For Tj^Q) , we have already chosen set of coordinates such 
that the induced metric on the string is conformally flat 
and X ■ X' — 0. In the same gauge, the solution to the 
Nambu-Goto equations of motion from eq. ([6|, [d^ — 
= 0, can always be expressed as a superposition 
of null left- and right-movers, 

X^(i,a) = ^(X^[a+]+X^[a_]) (11) 



and 



= ('T±,X±(a±)), X'l = 0- (12) 



where from henceforth we shall switch to using the vari- 
ables 



a± = t ± a 
and their associated derivatives 

d± = dlda±, 



(13) 



(14) 



instead of {t,a). Moreover it is advantageous to exploit 
the periodicity of the trajectory, i.e. X^^\t,a ± i] = 
X^^[t,(j\, to Fourier decompose the stress energy of the 
cosmic string, eq. (10 1, into the form 



2Mpi (NG) 



4Mpi ^ 



K) 



X (27r)(5 ka 



(15) 



(We denote U^^'V^ = U'^V + WYt".) The C^^^^f con- 
tain information about the trajectory of the cosmic string 
loop and are given by 



d(T-i 



exp 



d±Xl, 



where 



= p' - sp. 



(16) 



(17) 



For the scattering problem at hand, s — 1, so that the CiS 
become their complex conjugate under the interchange of 
p O p': C(j.\^)[p,p'] — C(-|-|f)[p',p]. (The overbar denotes 
complex conjugation.) For the pair production process 
described in |2], s = — 1. 

Let us denote the polarization vector of the incoming 
photon as [p] and that of the outgoing photon as e'^ [p'] . 
In the Lorenz gauge {d^^A^ = 0), which we shall use here. 



they satisfy e • p = e' • p' = 0. The stress-energy operator 
of the photons in eq. ^ , Wick-contracted with the state 
\p) representing one incoming photon of momentum p 
and with the state (p' | representing one outgoing photon 
of momentum p' is 



£ . + p . £{"g'/3} (18) 

-p{"e'«e-p'-p'{^e">e'-p 

aB I I -I -I I \ \ 

— Tj [p ■ p e ■ e ~ p ■ e p -e) I. 



The propagator of the graviton that follow from eq. 
(H) is 

F>aPfi^{k) = ^ iVafj^Vi^lS + Vai^VfJ^IB " VfJ^i'Vap) ■ (19) 

At this point, the quantum amplitude of a single photon 
with initial momentum p, scattering off the cosmic string, 
and acquiring a final momentum p' is then given by 



iMip^p') = -87:G^Tl^^^)D^,^f!T; 



(20) 



(Because the classical stress-energy tensor of photons 
is traceless in 4 dimensions, one would find that 

rj,fiT^^^/k\) This amph- 
tude has a diagrammatic representation shown in Fig. [T] 
We proceed to calculate |A^P, summed over all initial 
and final photon polarizations. As a consequence of Ward 
identities obeyed by electromagnetism, this amounts to 
the replacements ejje[,,e^,e^ —Vt^i^- From eq. ( [l5| , 
we will also face the product (27r)^(5(fco — (47r£/L))(5(fco — 
{4tt£'/L)) and their associated sums over all integers £ 
and £' . This product is zero unless the arguments of 
the two (5-functions are simultaneously equal to zero; this 
collapses the double sum into one, while the product of 
S functions become 



((27r)J(fco - {A7r£/L))y 

(27r)5(fco - (47r£/L)) (Total time elapsed). 

Define the amplitude-squared per unit time as 

\M{p 



\Aip ^ p': 



E 

polarizations 



p . 



where 



27r5(0) 



r+r/2 

27r(5(0) = fim / dt 

J-r/2 



(21) 



(22) 



/2 

Suppressing the £ index on the C(±|f)S, we arrive at 

4:TT£\ 



(23) 



\A\' 



OO 

E 



(24) 



with 



3 



ip' ■ vY 



- Hp' ■ p){{p' ■ c+W ■ cZ)(c^ ■ c_) + ip' ■ c^)ip' ■ • cl) 
- ip' ■ c^w ■ cZ){c+ ■ c_) - {p' ■ c+w ■ c^)(c^ ■ c^) 

+ \C^\'\p'-C+\' + \C+\'\p'-C^\'} 
+ 2{p'-pr{\C+-CZ\' + |C_nc+p - \C+ ■ |C±p ^ C^civap 



(25) 








FIG. 1; Feynman diagram for the amplitude of photon scat- 
tering off a cosmic string in motion. The ® represents the 
stress energy operator of a Nambu-Goto string in Eq. (10 1. 
The dashed hne is the graviton propagator in eq. (191. The 



rest of the diagram is the stress-energy of the photons in eq. 
( 18 1; the wavy hnes are the incoming photon with momentum 
p and the outgoing photon with momentum p' . (Drawn with 
JaxoDraw j6j.) 



where s is the same parameter occuring in (17 1. We 
display this dependence of \A\'^ on s so as to compare 
against the pair production result in eq. (26) of The 
primary differences going from the scattering to pair pro- 
duction amplitude are that, s changes sign from 1 to —1, 
and the k = p' —p in the C'(-|-|£) is replaced with k = p' +p. 
For the pair production case, since poiPo — 0, the sum 
over £ only runs over and the positive integers; whereas 
for the scattering case, i runs over all integers and zero. 
Whichever k = p' ± p is being used, k ■ C(^|-|-)(fc) — 
always holds; we have used this to arrive at eq. (25). 

With the scattering amplitude per unit time, we may 
insert it into the generic form of the Boltzmann equation 
for the scattering of bosons (see, for example. Chapter 5 

of H), 



dM^p'(27r)-= 



1 



d^p 



(2^)32po 



(26) 



(l^(p^p')l'/(i,p)(l + /(t,p')) 



W^P)! /(i,p')(l + /(^,P)) 



where d'*A'!^/dM'^p'(27r)^'^ is the number of photon scat- 
terings off the cosmic string in question, per unit time, 
per unit momentum space; and /(i,p) is the number of 
background photons per unit volume per unit momen- 
tum space. The first term describes the rate of photons 
entering the momentum bin p' due to their scattering 
from other momenta p; the second term is the rate of 
leaving the momentum bin p' due to the same scattering 
mechanism. 



One may check directly using eq. (25) that \Ai\ is in 



variant under the swap p -v^ p': p-p' is invariant, whereas 
as noted already, C±[p',p] = C±[p,p'], so p' ■ C± is re- 
placed with p-C± = —p' ■ C± (and p' ■ C± with — p' ■C±), 
where we have used k-C± = 0. Since every product of the 
p' ■ C±s and p' ■ C±s contains even powers of p', the minus 
sign can be dropped; hence, the swap p ^ p' amounts to 
the replacements (p' • C±,p' ■ C±) {p' • C±,p' ■ C±) 
in every term of The conclusion then follows 

if we observe that, the terms in {Ael"^ are either real 
or comes in pairs whose sum is real. This invariance 
|^^(p p')P — l-4f(p' ^ P)P means the quadratic-in- 



/ terms in eq. ( 26 ) cancel 



Furthermore, if we ignore the backreaction of this scat- 
tering process on the background distribution of photons, 
so that f{t, p) can be assumed to be time independent; 
then /(t, p) = 2(exp(|p|/r^) - 1)~\ the Bose-Einstein 
distribution for a massless spin-1 particle. For simplicity, 
let us also integrate over all possible final directions of 
the photon. Altogether, we are now lead to 



d'N^{t,\p'\) 
dM|p'|(27r)-i 

Ip I 



E 



2(27r) 



dipllpl 



(27r 



>3 



(27) 



x(27r)5 Ip'I-IpI- 



1 



— j IA(p^p)l 



1 



exp(|p|/T^) - 1 exp(|p'|/r^) - 1 



To be clear, our notation for the incoming p and outgoing 
p' momenta is 



P 



|p|(i,p), P' = Ip'I (i,??) 



(28) 
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so that, for instance d^p is the integral over the solid 
angle of the incoming photon. 

Note that, because scattering does not change the to- 
tal number of photons present, integrating the right hand 
side of eq. (27) with respect to |p'| (from to oo) would 
yield zero. We will later use this as a check of our nu- 
merical calculations. 

To summarize what we have found so far: for a given 
cosmic string loop trajectory (and the relevant left- and 
right-movers X±), first determine the modes C(-tms in 
eq. (16), and afterwards insert them into eq. (25) to 



obtainpl^p. Then the rate of photon scattering off the 
cosmic string per unit energy, is obtained by evaluating 
the Boltzmann integral in eq. (27). 



III. DEGENERATE KINKY LOOP 

In [2] , an analysis was performed for the gravitationally 
induced pair production rate of photons without special- 
izing to a particular cosmic string loop trajectory. This 
is possible because, for the pair production case, there 
exists asymptotic expansions of eq. ( 16 ) by utilizing the 
large dimensionless phase \k ■ X±\ ^ £ 1. (In other 
words, one can exploit the large ratio of the astrophysi- 
cal length of the string to the microscopic wavelength of 
the photon.) For instance, one may seek a saddle point 
of the C'±, i.e. search for configurations {p,p') where 
(p -I- p') ■ Xj. = 0. Because p,p' and X'^. are null vec- 
tors, that means p + p' must be spacelike or null for a 
saddle point to exist. That in turn implies p and p' are 
coUinear, since we require p ■ p' = |p||p'|(l — p- p') < 0. 
A slight deviation (in, say, the direction vectors p and 
pi) from this configuration would fail to satisfy the sta- 
tionary phase condition and therefore lead to a subdom- 
inant pair production rate; this in turn yields estimates 
on how very narrow regions (whose widths scale as some 
inverse power of £) in angular space provide the dominant 
contribution to the C± integrals. In all, what [2] found 
was that the energy emitted due to the pair of photons 
was largely independent of the mode number in other 
words, equal power in photon pairs were emitted at each 
harmonic |p| + |p'| — 'iiri/L. 

For the scattering case at hand, such an asymptotic 
analysis appears to be more difficult. For example, the 
stationary phase condition now becomes {p' —p)- X'^ = 0; 
but with the minus sign instead of a plus sign (as was the 
case for pair production) , the requirement that p' — p is 
spacelike or null yields no useful constraint on the config- 
uration (j),p'). Another important difference is that, for 
the pair production process, because |p'| + |p| — Ani/L, 
the rate at a given i, without performing the analogous 
infinite smu in eq. (27 1, is already a physical quantity 



of interest; it is the power generated at a given energy 
Airi/L. For the scattering process, on the other hand, we 
wish to ask how many net photons are being added or 
subtracted, at a given |p'|, from the background thermal 
bath. This means, for a fixed |p'| we have to perform 
both the solid angle integrals and the infinite sum over £ 
in eq. (27 1 in order to obtain a physical quantity. How- 



ever, at least for the two photon temperatures we have 
considered, we shall find that the dominant contribution 
to the infinite sum come from the small \£\ terms, specif- 
ically \£\ = 1,2, 3,..., 0(10- 100). 

So as to achieve a concrete estimate of scattering rates 
of photons off cosmic strings, we shall therefore, in this 
section, specialize to the "degenerate" kinky loop trajec- 
tory that was also used in T and [F. 

Specifically, the left- and right-movers of the kinky loop 
are 



X+(a+) 
X_(a_) 



(T+A < (7+ < f 



< L 



(7 B 



< cr_ < 



(L-cr_)B f < (T_ 



< L 



(29) 



where A and B are constant unit vectors. Now, denoting 



fcA = fcjA* 



we have, via eq. ( 16 1 



(/cA, -koA) 



C("_m-(-AB,fcoB) 



gi(feBi/8+7r£/2) 

L{kl - kl) 



8 2 



■ sm 



'-L+-£ 



(30) 



We first integrate over |p| in eq. (27), using the i5- 
function to set |p| = |p'| — Ati£/L in the integrand. Be- 
cause IpI > 0, we also need to introduce a step function 
0(|p'| — Airi/L) multiplying the result after integration. 
Next, for convenience, we set L = 1 and rescale 



(|p'|,|p|,r,)^47r(|p'|,|p|,r), 

followed by the definitions 

A^p'-p = p'-(|p'|-^)P- 



(31) 



(32) 



With the rescaling in eq. (31 ), the length |A|, as p-p' is 



varied from -1-1 to —1, ranges over 

|^| = ||p'|-|p||<|A|<|p'| + |p|=2|p'|-^- (33) 



Now, eq. (27) reads 



5 



256|p-||p|9(|p|) f 1 



\p\=\p'\~e 
sin^ (^^TT (^A • A 4 



er— 1 CT— 1 



(34) 



((A- A)2-^2^ 
(B- A)2-^2 



A'A^' (A'A^(1 + - 2£{\p'0'p - IpIp*^') + (5*^£^(1 - p • p')) 
gfcg/ (AfeA'(l + p . p') - 2£ (Ip'Ip'V - \P\P^P^) + S'^^i^i'^ -P-P')) 



Collinear limit The reader may worry that the pres- 
ence of {1 — p ■ j5')2 in the denominator of eq. (34 1 indi- 
cates the differential scattering rate may be singular in 
the collinear limit, a&p-p — > 1. This is not the case, but 
rather, 



p^p dM|p'|(27r)-id%/d% 
256|p'||p|e(|p| 



(35) 



(GnA*)^ 



_4/>4 I M -, Ip'l 

sin^ (2^ (1 _ A • p')) sin^ (2^ (1 - B • p')) 



The third line of eq. ( 35 1 contains all the remaining 



dependence on angles and is finite for all A, B, and , 
because I is an integer and sin(x)/x is finite at a; = 0. 

£ = Q Can Be Dropped We also note that the ^ = 



term in eq. (34) can in fact be dropped because, in this 
limit, the Bose-Einstein factors on the first line of eq. 
(34 1 goes to zero while the rest of the integrand is finite. 



Explicitly, the Bosc-Einstcin factors, when \£\/T <^ 1, 
can be expanded as 



lim 



e T — 1 




£ coth ( ^ 



^ \ 3 (cosh(l§ii 




2 f3csch2 



(36) 
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5 ) coth ( ) csch 



2 

2T 



48 



(For later use, we have expanded up to 0{£'^).) This leads 
us to the following expansion of the small i terms of eq. 



(Hi: 



_ (1^ 
2 \2) 



\l-p-p'J 7r4 T(elP'l/^-l 
/ sin (f Ip'IA ■ (j7 - p)) sin (f |p'|B ■ (j7 - p)) 



^\p'\A ■{p'-p) 
0{f). 



flp'IB.(p'-p) 



(37) 



For p ■ p' 1, this expression is non-singular for all |p'|, 
and therefore zero when the limit £ — > is taken. We also 
remark that, in the preceding discussion on the collinear 
limit, if we took p ■ p' ^ 1 followed by £ — 0, we would 



find the Boltzmann integrand in eq. ( 35 ) would read 



lim 



(Gnm) 



2 256|pf 



M 
£e T 



T(e- 







A 




(■k£^ 




B 


9 


it 




+ A 


9 






+ B 


Pl 



(38) 



(l + 0(£)) =0 



1 and I 







In other words, the two limits p ■ pl 
commute and yield a zero integrand. 

Average Over Loop Configurations If we were able 
to integrate over the solid angles with respect to the in- 
coming and outgoing photon directions (p and 1?), the re- 
sulting expression must only depend on A and B through 
the Euclidean dot product A-B, since the scattering rate 
is a scalar and there are no other vectors in the current 
problem. The result would allow one to see how the scat- 
tering rate per unit energy |p'| depends on the kinky loop 
configuration. However, in this paper, we shall be con- 
tent with obtaining an estimate of the rate of scattering 
averaged over all possible kinky loop configurations - if 
we assume a uniform distribution, this means we shall 
now examine 

j2 



/ d^A^,(t,|p'|) 
\dM|p'|(27r)-i 



§2 



df^A 

47r 



§2 



d»B d"jV^(t,|p'|) 

An dM|p'|(27r)-i' 
(39) 



Notice the Boltzmann integrand in eq. (34 1 has factor- 
ized into three portions; the first line is independent of 



6 



A and B; the second depends solely on A and the third 
solely on B. This means the averaging in eq. (39) also 
factorizes into two independent solid angle integrals, one 
with respect to A and the other with respect to B. In 
fact, the primary integrals wc need are 



f3 



cm, 

g2 47r 



^sin2(i^( 



A- A- 



(A • A)2 - £2 



An 



B'B^ 



(40) 



(B- A)2-^2 



By a flip of parity, A -> — A or B ~> — B, we see that 
these two integrals are actually the same, rj = Ig = 
Therefore let us focus on the first. Observe that the 
only vector in the integrand is A, and the integral trans- 
forms as a rank 2 tensor under 0(3) transformations: for 
A^ -> U'jA^, where U e 0(3), we have U'jU\PK 
Therefore we may fix the form of P^ to be 



5'^ A, 



(41) 



where Aa, As are 0(3) scalars. Contracting both sides 
with 5ij and with A^A^ allow us to solve for and A^ 
in terms of the two scalar integrals 



In = 



'A- A 



47r 



(< 



A- A) 



<'2 



, (42) 



+1 sin^ (i7r(|A|c + ^)) 
dc- 



(A2c2 - £2 j 



where c = A • A/|A|. The results are 

Aa = 
As = 



3/i + 3f2 /2 - A2 /2 

2A4 

A2 h-h- P h 



2A2 



(43) 
(44) 



and 



2/i 



2^|A| 







l^ln 


£+|A| 



Ci(^|^-|A||)-Ci(7r|£+|A| 



2h 



4^3|A| 



In 



^+|A| 



lAI 



(45) 



+ 



(Ci(^|^+|A||)-Ci(7r|£- |z 



cos(7r(|A| +£)) - 1 

2^2(£2 _ A2) 



The Si and Ci are the Sine and Cosine Integral respec- 
tively. Whereas the In minus Ci terms may also be re- 
expressed in terms of Cin and the Euler-Mascheroni con- 
stant 7e (see [5]), 



Cin(z) = ln(z) - Ci(z) -I- 7e, 



(46) 



which is an even power series - this is why we may put 
an absolute value around the arguments of In and Ci. 



With the Ii,2 in eq. (45), and recalhng |p| = |p 
A = p' — (Ip'I — l)p, the Boltzmann integral in eq. (34) 
averaged over loop configurations now becomes 



^^N,{t,\Yy'\) \ , ^ 512|p'||p|e(|p|) f 1 

dMip'i(2.)-r(^-^) 

IpNIp'I- 



1 7r2A4(l-p.p')M -1 e^-1 



(47) 



A* A J 



J2 (3/1 + I2 - A2 + (a2 h-h- e h) 5'^ 
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where c' = p - f/ . Even though p and f/ appear in the in- 
tegrand, after all the Euclidean dot products are carried 
out, the integrand in eq. (|47]) will only depend on them 
via the combination p ■ p' . Again using 0(3) invariance, 
we have reduced the Jga dfip integration to 2tt times a 
one dimensional integral with respect to c', after the re- 
placement p ■ — > c'. After this integral with respect 
to c' there will not be any dependence on p nor p', and 
therefore we have replaced the solid angle integral with 
respect to jf with in. 

Zero Ip'I limit Taking the zero outgoing energy 
limit, Ip'I 0, of Eq. pTl) we get 



dM|p'|(27r)-i 



6ALT^{GN^l)^ ^ Cin2(27r^) 



-220(Lr^)(GNM) 



2 



£3 



(48) 



As we shall soon see, the rate of scattering (away from 
Ip'I = 0) will continue to be negative for low energy 
photons, before becoming positive at higher energies. 
The motion of cosmic strings kick photons from lower 
to higher energies. 



We will now evaluate eq. ( 47 ) 
10^2 and 10^6. For a photon 



Change-of-variables 
numerically for T — 
temperature of 3 K, this corresponds to loop sizes of 
L{T = 10^2) = 3 kpc and L{T = 10^6) = 3 x 10* kpc 
respectively. (For comparison, the diameter of the Milky 
Way is 30 kpc.) Let us change our variables from c' to 
A = |A|, i.e. 



where — AnT in this formula denotes the original 
(un-rescaled) photon temperature (see Eq. (31)), and we 
have reinstated L. This non-zero scattering rate indi- 
cates that, with a sufficiently high cosmic string density, 
there could be a significant backreaction on the low en- 
ergy photon distribution. 

A| 



To get to eq. ( 48 1 we have set 



In particular, starting from eq. (45) 



-> 1^1, and 
as Ip'I -> 



h 
h 



Cin(27rf) 
Cin(27r^) 



- Til Si(27r^) 



8^4 



(49) 
(50) 



,'|2 



|p'|-£)2-A2 



dc' 



to obtain 



2|p'|(|p'|-^) 
AdA 



Ip'KIp'I-^)' 



(51) 
(52) 




d^iV,(t,|p'|) \ .2 V 89(|p'|-£) / 1 I 

dtd|p'|(2^)-i/ ^ l^e^-^-l e^-iy^m 

X [ J (Gin (^(A - t)) - Gin (7r(A + t))) (SA^ - A^ (^^ _,_ (2|p'| _ £)2) + U''{2\p'\ - if) 



(53) 



(3(2|p'| - If - A^) (2A ((-1)^ cos(7rA) - l) + 7r(A2 - f) (Si (7r(A -I- (.)) + Si (7r(A - I)))) 



Let us understand qualitatively the integrand in eq. ( 53 ) , 
over the range of integration, |^| < A < 2|p'|— £. GiriyKz) 
is an even function of z, starting from zero at z = and 
rising, for large |z| ^ 1, to go as ln(7rz) (see eq. ([46|). 
Gonsequently, the l'^{Gin (7r(A - £)) - Gin (7r(A -hI))} 



in eq.(53) is an even function of peaks at A = \l\ and 
proceeds to fall off as (2/A)(cos(7rA) - 1) + 0{{l/J\f) 
for A ^ Hence, the second line of eq. (53) grows at 

Whereas 




most as a cubic polynomial in A for A 3> 
Si(7rz) is an odd function that rises from zero at z 
and asymptotes to 7r/2 for large ttz 2> 1; and thus the 



Si(7r(A + <')) + Si(7r(A-£)) term asymptotes to tt for 
A £. We see that the third line of eq. ( 53 ) (inside the 
(...)2), grows as a quadratic polynomial multiplied by 
3(2|p'| -£)2 - A^, which decreases as A ^ 2|p'| from 



below. With the 1/A^ (on the first line in eq. (53)) mul- 
tiplying the entire squared quantity, we therefore expect 
the integrand to peak close to the lower end of the inte- 
gration A ^ |£| and fall off with increasing A. Moreover, 
for fixed |p'|, large |^| ^ 1 terms ought to be suppressed 
relative to their smaller |£| counterparts. These observa- 
tions will aid us in the ensuing numerical integration - in 
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FIG. 2; The number of thermal photons scattering off the 
kinky cosmic string loop (of size L) described in eq. (29 1, per 



unit time, per unit energy. There are actually two curves here, 
one at temperature 10^^(47r/L) and the other at 10^^(47r/L); 
because for |p'| ~ the rate scales as (LT-y)^, they coincide. 

Appendix we discuss some of the technical issues we 
face. 

We remind the reader that the rescahng in eq. (31) is 



in force from equations ( 32 ) through ( 53 1 , except for eq. 

In Fig. [2] we present the results for eq. (53), when 
= (47r/L)1022 and {4tt / 1)10'^'^ , with respective maxi- 
mum scattering rates per energy (roughly) ±10'*^(G'N/i)^ 
and ±10^^(GnA*)^- That the higher temperature yields a 
larger rate is what one would expect on physical grounds, 
since there are more photons available per unit volume 
at higher temperatures. (Photon number density is pro- 
portional to the cube of its temperature.) Since temper- 
ature enters the current problem only through the Bose- 
Einstein factors in eq. ([27]), and since it is these Bose- 
Einstein terms that are responsible for the existence of 
positive and negative scattering regions, we see that the 
energy scale over which this transition from negative to 
positive rate takes place is itself of the 0{T^). 



IV. SUMMARY AND DISCUSSION 

We have calculated the amplitude for inelastic scatter- 
ing of photons off a cosmic string loop. We have used the 
result to calculate the spectral distortions from a loop of 
string placed in a thermal background of photons. The 
shape of the spectral distortion is shown in Fig. [2] 

It is of interest to estimate the fraction of photons that 
will be inelastically scattered by a distribution of cosmic 
string loops. First we note that the scattering rate shown 
in Fig. [2j can be estimated as 



dtduj 



(Gnm)'(LT^)', uj^T., 



(54) 



the frequency of the photon. Moreover, here we are only 
concerned with the scattering when a; ~ T-y because this 
spectral region contains the most number of photons. 
The important scaling as has been verified numeri- 
cally by comparing the two plots in Fig. [2j it can also be 
checked by using rescaling arguments on Eq. ( 34 1 . 



Now we can multiply the scattering rate by the loop 
number density to get the number of scattered photons 
per unit time per unit volume. Since the scattering rate is 
proportional to L^, the largest loops will be most impor- 
tant for this process. Such loops have size comparable to 
the horizon and so we take L ^ t, whereas the loop num- 
ber density is approximately where t is the cosmic 
time. Then we can integrate over time to get the total 
change in number density of photons Art per frequency 
interval, 



d(An) 
duj 



tinitial 



i3 



"'tinitial ^ 



(55) 



where Ti is the initial temperature from which we start 
counting the number of scattered photons. (We have used 
the fact that the photon temperature scales inversely as 
the cosmological scale factor, T-y c>c l/a(i), assumed a{t) 
goes as some power law, and the initial scale factor is 
much smaller than its present value a(iinitiai) ^ a(inow)-) 
To get the fractional change in photon number density 
at some frequency, (d(An)/dw)/(dn/daj) ~ An/n evalu- 
ated at the present epoch, due to scattering by strings, 
we divide Eq. ( 55 ) by the number density of ambient 
photons at the present time with frequency w ~ Tp, i.e. 
dn/doj 



T§. Then 



An 



2 „2 



{G^f^rz 



(56) 



where in this section we are using uj — |p'| to denote 



where z is the redshift after which the scattered photons 
do not get re-thermalized by Compton scattering. We 
will not go into a detailed calculation here but simply 
mention that z is definitely higher than the recombina- 
tion redshift (- 10^) and possibly ~ 10'' [3]. 

With GnM ^ 10^'' and using z « lO'', the spectral dis- 
tortion from inelastic scattering of photons by the gravi- 
tational field of strings is at a level ^ 10~^. The distor- 
tions are similar to the ?/— distortions of the CMB in that 
the total number of photons is unchanged. However, the 
shape is different and almost completely resides in the 
low frequency part of the spectrum. 

As in the case of gravitational wave emission jlOj . 
cosmic strings produce a burst signal in addition to a 
stochastic signal. It is possible that cusps and kinks on 
strings can up-scatter ambient photons in a burst, and 
that this may provide another observational signature. 
We leave the exploration of this possible signature for 
future work. 
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Appendix A: Evaluating the Boltzmann Integral 
Numerically 

In this appendix we describe some of the technicalities 
involved in the numerical evaluation of the Boltzmann 
integral in eq. ( [53| . To begin, let us re-express eq. ( 53 ) 
in the following schematic form: 



± (Al) 



dM|p'|(27r)-i/ ■ ■ 

2\p'\~-e 



with 



1 



1 



(A2) 



One may attempt to use the NIntegrate function in Math- 
ematica to calculate numerically the dA-intcgral portion 
of eq. (All, for a given |p'|, over as many £s as possi- 
ble, before summing them up. One would find that the 
l-^l ~ C(l — few) terms are the dominant ones, but these 
same terms are such that, the pair zL\£\ yields dA inte- 
grals of nearly equal magnitude but opposite signs. It is 
therefore prudent to combine the positive and negative £ 
terms, for \£\ < |p'|, before integrating, i.e. 



d^N^it, Ip'I 



dM|p'|(27r)-i 



+ 00 



(2|^)J 



e=i 



fcLlp'IJ+i 



where 



S(2|£) 
S(3|£) 



/ {Efle + S-e^-i) dA 

Ji 

■2\p'\+t 



2\p'\-e 

2|p'|+^ 



(A3) 
(A4) 
(A5) 



Notice we have also changed the summation variable 
such that now £ is summed over all non-zero positive 
integers. 

For , it is often useful to expand the Bose-Einstein 
factors using the small £ limit in eq. ( 36 ) before applying 



NIntegrate. Specifically, by decomposing the integrand 
into ^-even plus ^-odd terms, and denoting X± = ±I_|, 
we have 



= |P'I/T 



(IK) 



T (elP'l/T^ - 1) 



2\p'\-i 



dA 



(A6) 



coth 



{^2t) 



T 



-I- 



2 (3csch2 (^) 



cosh 



12 
iP^i 

T 



+ 5 coth 



iP^ 

2T 



csch 



lp!i 

2T 



48 



Of course, there are only two distinct integrals here. 



namely 



2|p'h 



I±dA; the temperature dependent pieces 



of do not take part in the integration. What we are 
emphasizing is the cancellations that occur between the 
£-odd and £-even portions of Te when summing over the 
originally positive and negative £ terms, at each order in 
the £/T expansion. For all the cases where we used this 
expansion in cq. ( A6 ) , we have found that the sum of the 



third and fourth lines is always much less than the second 
line - usually much less than a percent - and therefore 
may be neglected. 



each 



> 10, we evaluate B, 



(1,2K) 



for £ 



For 

1, 2, 3, . . . , 10 and sum them up. For £ > 10, we per- 
form the integrals for ten £s per decade, for example, 
^ = 20, 30, . . . , 100 and £ = 2000, 3000, . . . , 10000, up to 
a few orders of magnitude beyond that of T. (For small 
Ip'I ^ T, due to numerical underflow we do not inte- 
grate beyond £ ~ T but we simply stop at some £ that 
yields B(^i 2^3\e) that are many tens of orders of magni- 



tude below B, 



(1,2|£=1)- 



Then we use the Mathematica 



function Interpolation to compute a best fit line to the re- 
sulting numerical data and proceed to approximate the 
summation for £ > 10 by an integral over the fit range, 

J2i=Li ~^ //=io^^- ^^^'^ that, for each |p'|, the sum 
over the first ten £ = 1, 2, . . . , 10 is usually responsible 
for roug hly 91 percent of d^7V^(|p'|)/dM|p'|. 

When numerically evaluating 6(1^2. 3|£) we have found 
it useful to not do it over the entire range of A, but 
rather start from integrating from the lower limit Ai to 
min[10Ai, A2], where A2 is the upper limit (for instance, 
Al = £ and A2 = 2|p'| - £ for -B(i|^)), then from Ai 
to min[100Ai, A2], etc., until the error incurred is less 
than one percent. This algorithm is justified because, as 
already explained after eq. ( 53 ) , we expect our integrand 
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to always peak near A ~ \£\. Practically speaking, this 
saves Mathematica considerable amount of time and effort 
when sampling the integrand. 

As a check of our calculation, we fit a line through the 
numerical data points of our plots in Fig. [2] and find that 
the area above the |p'|-axis is equal to the area below 



the |p'|-axis at the sub-percent level. This is both con- 
sistent with our demand for percent-level accuracy in the 
numerical integrations of S(i 2,3|£), as well as the physical 
requirement that scattering processes do not change the 
net photon number. 
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